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We investigate the two-point correlations in the band spectra of periodic systems that exhibit 
chaotic diffusion in the classical limit, in terms of form factors with the winding number as a spatial 
argument. For times below the Heisenberg time, they contain the full space-time dependence of the 
classical propagator. They approach constant asymptotes via a regime, reflecting quantal ballistic 
nrotion, where they decay by a factor proportional to the number of unit cells. We derive a universal 
scaling function for the long-time behaviour. In the limit of long chains, our results are consistent 
with expressions obtained by field-theoretical methods. They are substantiated by numerical studies 
of the kicked rotor and a billiard chain. 
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I. INTRODUCTION 

A large class of systems, among them most promi- 
nently solid-state systems, are organized as repetitions 
of identical or near-identical units. If in such a system 
the classical dynamics is chaotic and the unit cells are 
connected, then globally this leads to a diffusive spread- 
ing of trajectories, irrespective of the presence or absence 
of static disorder. Regular motion on tori does not con- 
tribute to diffusive spreading: Closed tori do not allow 
for transport along the lattice, while open tori give rise 
to ballistic spreading. The spectral and transport prop- 
erties of extended quantum systems, in contrast, depend 
sensitively on the degree of translational symmetry. On 
short time scales, however, where the quantum dynamics 
still closely follows the classical, both periodic and dis- 
ordered quantum systems exhibit (apart from an initial 
ballistic phase which is of no interest for the following) a 
spreading of wavepackets with the characteristics of the 
chaotic diffusion in their classical counterparts. The sig- 
nature of this phase in the discrete spectra of disordered 
systems with Anderson-localized eigenstates has been in- 
vestigated previously [l]|J] , using techniques in the spirit 
of Berry's semiclassical derivation |TJ,H of the spectral 
form factors for the canonical random-matrix ensembles. 

In the present work, we report on a study of periodic 
systems with band spectra and eigenstates of Bloch form, 
focusing on the signature of spatial order and dynamical 
disorder in the spectral two-point correlations. We ex- 
ploit the existence of a second conserved quantum num- 
ber besides the energy, the quasimomentum, to define 
form factors in the canonically conjugate space spanned 
by time and winding number. They are related to the 
spatially coarse-grained propagator and therefore ideally 
suited to extract dynamical information from the band 
structure, without recurring to a local spectrum. Here, 
we evaluate this relation over the entire time evolution 
of the form factors. In the semiclassical regime, we show 



that the full space and time dependence of the classi- 
cal propagator is contained in the form factors. Chaotic 
diffusion in spatially periodic quantum systems can thus 
directly be identified in the band structure. Similarly, 
the spectral signature of the crossover to quantum ballis- 
tic motion can be analyzed on basis of the full quantum 
propagator. In the limit of a large number of unit cells, 
the form factors exhibit a marked peak in the vicinity of 
the Heisenberg time, a consequence of the clustering of 
levels in quasicontinuous bands. We obtain a universal 
scaling function for the long-time behaviour. 

Form factors containing the correlations of levels across 
the Brillouin zone have previously been studied using the 
supermatrix nonlinear sigma model flqjq] . Our approach 
is complementary in that it emphasizes the important 
concept of winding numbers, providing direct access to 
spatial information. In combination with semiclassical 
techniques, this allows to draw a particularly transpar- 
ent picture of the physics over all time regimes. We shall 
demonstrate that in the case of diffusive spreading and 
in the limit of an infinite number of unit cells, our theory 
allows to reproduce the principal results of the sigma- 
model approach. At the same time, however, we are able 
to go beyond those results in that we need not make any 
assumptions as to the number of unit cells in the system, 
nor on the mode of density relaxation. 

Our theory therefore applies to a broad class of sys- 
tems. They include semiconductor superlattices support- 
ing chaotic electron motion ("antidot arrays") [R-E0[, 
quantum-optical systems involving periodically modu- 
lated standing- wave fields pi] , and KAM systems with 
toroidal chaotic layers in phase space containing chains 
of regular resonance islands, as they occur frequently in 
molecules |1SJ| . For this type of systems, there exists a 
large body of spectral data, both experimental and nu- 
merical ]7|~[l2|]. By analyzing these data as explained in 
the following, the dynamical information encoded in the 
respective band structures can be extracted. 



II. THE GENERALIZED FORM FACTOR 

The setup we have in mind is a finite chain of N iden- 
tical unit cells, with cyclic boundary conditions at its 
ends. Spatial disorder within the unit cell is not re- 
quired. We restrict ourselves to quasi-one-dimcnsional 
lattices, since the generalization to higher dimensions is 
straightforward. 

As a consequence of periodicity, the spectrum can be 
decomposed into N subsets, each of which corresponds to 
one of the N irreducible representations of the group of 
lattice translations T(nae x ), n integer. Together, these 
subspectra coalesce into discretized bands and approach 
continuous bands in the limit N — > oo. Dynamical and 
spectral quantities specific to one of the irreducible rep- 
resentations M,H are constructed using the projectors 



N-l 



Pm = jr y, Xm(n)T t (nae a 



(1) 



onto the corresponding subspaces. They invoke the 
group characters Xm(n) = exp(in# m ). We refer to the 
9 m = 2nm/N, m = 0, . . . , N — 1, as Bloch phases. The 
symmetry-projected Green function is defined as 



G m {E) = P m G(E), 



(2) 



where G(E) is the Green function for the full chain. From 
G m (E), other Bloch-phase-specific quantities can be de- 
rived as if they pertained to the full spectrum of a sys- 
tem without spatial symmetry. For example, the Bloch- 
phase-specific spectral density is related to the trace of 
the corresponding Green function, Eq. (g), in the usual 
way (see, e.g., Ref. @), 

d m (E) =y2KE-E a (6 m )) = --Imtr [G m (E)}. (3) 



Here, the trace extends only over a single unit cell. 

The basic energy scale in the following is the inverse 
mean spectral density per unit cell, l/(d uc ), or equiva- 
lently, the mean separation of neighbouring bands. We 
define also the Heisenberg time with respect to the unit 
cell, £jj = 27r?i(d U c). Accordingly, we scale time as 
t = £/<h an d energy as r = (d uc )E. In these units, 
the size of the spectral window considered is Ar, roughly 
the total number of bands. 

The time-domain counterpart of d m {E) is the ampli- 
tude 



a m {r) 



dr e 



r d m (r/(d uc )). 



(4) 



By performing another, now discrete, Fourier transform 



with respect to 6 m []L4 17 1, which amounts to going from 
the Bloch-phase to the winding-number representation, 
we define the amplitude 



1 N_1 



» 



— ^ ex P(i™#m)amC 
m=0 

/ dq (q + nae x \ U(rtn ) | q) 

-/unit cell 



(5) 



where ae x generates the lattice, je^l = 1. Winding- 
number-specific form factors are defined as 



K n ( T ) = -L\~a n (r)f 



(6) 



Substituting Eq. (pj) shows that the K n (r) comprise pairs 
of levels with all Bloch phases. In particular, K (t) cor- 
responds to the form factor for the entire spectrum, irre- 
spective of spatial periodicity. Eq. (0) represents a partial 
trace of the propagator U. The K n {r) can therefore be 
interpreted as probabilities to return after encircling the 
unit cell n times. 



III. SEMICLASSICAL REGIME 

A semiclassical account of the symmetry-projected 
spectral quantities is achieved on the basis of a gener- 
alized concept of periodic orbits |13 lljJlT]]. It becomes 
transparent in a symmetry-reduced representation of the 
chain, where the two boundaries of the unit cell con- 
nected by the translational symmetry are identified. In 
this way, the unit cell assumes the topology of an annulus, 
possibly times additional dimensions. An orbit periodic 
in this reduced space can be classified according to its 
topology, expressed by its winding number, the number 
of times it runs around the unit cell before closing. The 
contribution of a periodic orbit j to d m (E) contains the 
additional phase —n^9 m , where rij is the winding num- 
ber of this orbit [13[l^,|lj. By the two Fourier transforms 
that lead from d m (E) to d n (T), we obtain the semiclas- 
sical trace formula 



hi(r) = Y^ 



rf (.S 3 {E) 

, exp i J I 

, \/|det(M,-l)| V n 

X 8l/Ar{T — Tj)^(n-n 3 )modAT, 



^2 



(7) 



(p) 



with t vp ' 1 , Mj, Sj, and fj,j denoting, respectively, primi- 
tive period, monodromy matrix, action, and Maslov in- 
dex of the periodic orbit j. Besides the usual amplitude 
and phase factors, this trace formula has attained two 
delta functions: A broadened <5i/Ar( r— T j) of width 1/Ar 
picks out orbits with scaled period Tj ps r. An A-pcriodic 



Kronecker delta S 



(n — rij ) mod AT 



selects periodic orbits with 



a winding number that differs from n at most by an in- 
teger multiple of N. 

According to the above interpretation, the K n (r), for 
t < 1, should be related to the classical probabilities 
Pn (t) to return, after n windings, in the symmetry- 
reduced phase spac e. Indeed, within the diagonal ap- 
proximation plpljl^], we derive 



X(-)(r)= 7 „rP( cl )(ri H ), 



T< 1. 



(8) introducing the dimensionless parameter g u 



N 2 



Here, we have neglected the contribution of repetitions 
of shorter periodic orbits. We have not taken the oc- 
currence of self-retracing orbits into account in order to 
replace individual degeneracy factors, expressing time- 
reversal (T) invariance, by a global 7„. Reflecting weak 
localization as a function of n, it takes the value 2 if orbits 



with i 



n are generically T degenerate and 1 otherwise. 



The winding-number representation thus enables a direct 
and natural access to weak-localization enhancements in 
the form factor. In the Bloch-phase representation, by 
contrast, weak localization is reflected in a smooth tran- 
sition from GOE statistics near the symmetry points of 
the Brillouin zone to GUE statistics elsewhere 0H0,[18|. 
In order that an orbit contribute to Pn (t) , it must 
be periodic up to a lattice translation by nae x . Assum- 
ing that the long periodic orbits spread as the generic, 
non-periodic ones, we express Pn (t) in terms of the full 
classical propagator p(r', r; t) as a partial trace, 



Pi cl) {t)= I drp(r + na(0,e x ),r;t) 

-/unit cell 



(9) 



where r = (p, q) denotes a point within the unit cell on 
the energy shell H . Equations (ph , (Bh show that the gen- 
eralized form factors, in the semiclassical regime, relate 
the band structure to the full, if coarse grained, classical 
propagator. 

The validity of Eq. (pft is not restricted to any spe- 
cific form of relaxation of the classical distribution, pro- 
vided the underlying classical dynamics is predominantly 
chaotic. For example, billiard chains connected only by 
narrow bottlenecks show a marked deviation from nor- 
mal diffusion on the time scale of the escape from a single 
cell. The generalization to higher-dimensional lattices is 
straightforward. Also there, P n ° (t) takes forms signifi- 
cantly different from diffusion in one dimension. 

As a specific example, we evaluate Eq. (pf) for nor- 
mal diffusion in one extended dimension. For an in- 
terval of length L = Na with cyclic boundary condi- 
tions, the diffusion equation is solved by the propagator 
p(x\x;t) = £( modL )(x'-:r,Z^),where£( modp )(x,cr 2 ) de- 
notes a normalized Gaussian of period p and variance a 2 . 
We assume that in the nonperiodic dimensions of phase 
space, the relaxation towards equidistribution is rapid 
on the relevant time scales. For times t <C id = L 2 /irD, 
the Thouless time for the full chain of length L, diffu- 
sion is free, p(x', x; i) = {2irDt)- 1 / 2 exp(-(x'-x) 2 /2Dt), 
while for t ^> id, equidistribution p(x' , x; t) = 1/L is ap- 
proached. If, as it is the case here, p(x! , x; i) = p(x' — x; i), 
the partial trace of the propagator amounts to multipli- 
cation by the cell size. We find for r < 1, 



n V ' N \N 



= / ln\/T/2g u 
InT/N 



yuc' 

irN 2 

7rn 2 /2g 1; 



r « N 2 /g u 
t » N 2 /g u 



tn/td = 2n 2 Ti(d nc )D/a 2 . Since we do not require static 
disorder and diffusion within the unit cell, the interpre- 
tation of <? uc as a conductance is purely formal. 

With respect to g nc and N, we distinguish two regimes: 
For g uc 3> N 2 , the classical dynamics becomes ergodic 
before the energy-time uncertainty relation allows to re- 
solve the inter-band spacing l/(<i uc ). The sampling by 
the discrete Bloch phases is then too coarse to reveal the 
continuous bands underlying the discrete levels, and the 
full spectrum appears as a superposition of N indepen- 
dent spectra. The condition g uc ^> N 2 can also be ex- 
pressed as 2n/N 3> CO rr, where 2n/N is the Bloch-phase 
spacing, and 6 COII = 2ir ^J 'ir / g uc is the spectral correla- 
tion length with respect to variation of the Bloch phase 
|| . Only for g uc <C A^ 2 , the arrangement of levels in 
bands is felt by the two-point correlations. A number of 
N9 cott /2tt levels then contribute coherently to the spec- 
tral correlations on time scales r < 1. In this case the 
second option in Eq. din) , the ergodic regime of the clas- 
sical dynamics, is irrelevant. In the space-time domain, 
this amounts to the diffusion cloud still being well lo- 
calized within the chain at r = 1. Increasing the chain 
length beyond N = ^/g uc merely results in a finer resolu- 
tion of the bands. 



IV. QUANTUM BALLISTIC REGIME 

Equation ( |f0| ) was derived using the diagonal approx- 
imation with respect to the classical phases. The peri- 
odic orbits occurring in the underlying trace formula are 
those of the symmetry-reduced space. The break time 
beyond which the diagonal approximation ceases to be 
valid is therefore the Heisenberg time for the unit cell 
£h, or equivalently, r = 1. This means that Eq. (10) de- 
scribes only the spectral correlations on scales of a typical 
inter-band spacing or larger. Therefore, we adopt a dif- 
ferent approach to K n (r) for r > 1, corresponding to 
energy scales of the inter-band spacing and below. Start- 
ing anew from the definition (0), we use Poisson resum- 
mation to replace the sum over m by an integral over 



a n (r) 



— V / d e Kn^N)ey e -2nir a (9)r (n) 



(10) 



where r a (8) — (d uc )E a (9). For large r the phase of the 
integrand is rapidly oscillating, and the integration can 
be performed within an asymptotic approximation. Pro- 
vided the bandwidth is of the order of the inter-band 
spacing, this approach is justified for times r ~ 1. It 
will in fact be seen in Section M that this approximation 
exactly reproduces the large r behaviour in the limit of 
N — > co. For simplicity we disregard special cases like 
inflection points or higher-order extrema which can be 



treated, e.g., by a uniform Bessel- function approximation 
JL9| . Saddle-point integration leads to the condition 



2*r , a (0 j (y)) = v 



(12) 



for points 9j(v) of stationary phase, with v — (n — 
fj,N)/r. Replacing dimcnsionless quantities by unsealed 
ones gives the equivalent condition v a (k) = (n — fiN)a/t, 
where v a (k) is the group velocity for the band a at quasi- 
momentum k = h9/a. It expresses the ballistic motion 
of Bloch waves. For r > 1, the phases 2nir a (0j)T left by 
the saddle-point integration can be considered random. 
Upon squaring to obtain the form factors, we therefore 
drop the off-diagonal contributions and get 



oc 

K n (r) = 1 -Y. 

/J,— — oo 



F 



n - fiN 



(13) 



where F[y) = (47r 2 Ar)- 1 E a E i KC^M)! -1 is a P os " 
itive function depending on n and r only through (n — 
IxN)/t. Similar to what we found for r < 1, the K n (r) 
can be interpreted as sums of winding-number distribu- 
tions F(i/(n, t)), each of which is stretching from a center 
at n = fiN, but now as a linear function of time, so that 
the variance a 2 grows quadratically. They are normalized 
by the common prefactor 1/r. In order to determine the 
as yet unknown function F(iA_, we use a heuristic argu- 




ment in the spirit of Refs 
the semiclassical expression 
t = 1. Expanding 



" k ' N y KN'ttN 2 ) 



and extrapolate both 
and Eq. ( |l3|) towards 



E expf- ^:^" ) 1 141 



2ffuc ^ V 2r 9 

/i= — OO V 



we match Eq. @ with Eq. @ at r = 1, 



OC 1 OC 

E F(n - "*> = -755= E 



exp 



/l— — oc 



7r(n - /iiV) 2 



2ff u 



(15) 



Comparing both sums term by term, we obtain F(y) = 
{2g nc )~ 1 / 2 exp(— ttv 2 /2g uc ) and thus, for the regime r > 
1, 



K n {r) 



7" C (n 

N 



odl) 



n g u 



N' irN 2 



(16) 



In the presence of additional symmetries, the spectral 
statistics can be that of the GOE throughout the Bril- 
louin zone, possibly with weak localization enhancements 
in the vicinity of the symmetry points (cf. the discussion 
of the quantum kicked rotor below). In this case, an 
analogous matching procedure applies. 

For the long-time behaviour of the entire set of K n (r) , 
Eq. ( fl6| ) implies the following scenario: Initially, Ko{ T ) 



decays as 1/r as long as only the term with \x = con- 
tributes significantly. As terms with larger \x attain a 
comparable magnitude, all K n (r) approach an asymp- 
totic constant "/ n /N (an exact evaluation pl| gives a 
correction —2/N 2 in the presence of the band symme- 
try d- m (E) = d m (E), if N is even). The asymptotic 
domain is reached atr« N \J 'n / g uc - This is the effective 
Thouless time for ballistic spreading. It corresponds to 
the time when the uncertainty relation allows the resolu- 
tion of the typical separation 2ir/(N6 corY (d uc )) of neigh- 
bouring discrete levels. It will be shown in the sequel 
that the expression ( |l6| ) reproduces the numerical data 
surprisingly well for r > 1. 




FIG. 1. Space-time dependence of the theoretical pre- 
diction for K n (r) according to Eqs. dlG) and ( |16| ) for 
AT = 512. 

Fig. u\ summarizes our results for the time depen- 
dence of the K n {r). We can identify three distinct time 
regimes: a semiclassical phase where the form factors 
reflect the classical diffusive spreading, a regime dom- 
inated by quantal ballistic motion, and an asymptotic 
regime reflecting the discreteness of the spectrum on the 
finest energy scales. These long-time asymptotes vanish 
for large N as 1/N. At the same time, the maximum of 
Kq(t) near r = 1, all other parameters being kept fixed, 
is proportional to 9 C on and independent of N. Therefore, 
in the regime g uc ^C iV 2 where proper bands exist, there 
is a crossover during which Kq(t) decays by a factor of 
the order of N9 COTI . Thus, together with its rise in the 
initial, semiclassical regime, Kq(t) attains a peak in the 
vicinity of r = 1 which expresses the clustering of lev- 
els into bands 0. The stationary-phase condition JI4), 
with fi = (only this term contributes near r = 1), shows 
that the peak is associated with the extremal points in 
the bands, that is, with the van-Hove singularities |2(J in 
the full spectral density d(E). 

Finally we note that for small N, in particular for 
N = 2, the coarse-grained density relaxation does not fol- 
low a diffusion law if there is no static disorder within the 
cells. An evaluation of the semiclassical and the quantum 



domains along similar lines as sketched here, gives access 
to the statistics of tunnel splittings in double billiards. 



V. THE LARGE-A LIMIT 

In the present section we discuss the limit of N — > oo. 
In this limit, our results simplify considerably, since only 
the term /i = contributes in Eqs. ( E3 ) and (O) ■ For 

short times r ~ 1 the K n {r) can be replaced by a single, 
diffusively spreading Gaussian as given in the second line 

of Eq. dlQ). In the long-time regime r ~ 1, we obtain 



K n {r) = 



In 



VWu 



exp 



2g u 



(17) 



This implies in particular that Kq(t) = l/y / 2g uc T. 

The same limiting behaviour has been calculated 
within the framework of a nonlinear sigma model pi . We 
shall now compare Eq. ( J17| ) to the corresponding results 
in Ref. pj . We make use of the fact that the Bloch phase 9 
can be viewed as a generalized Aharonov-Bohm flux and 
in the present limit, becomes a continuous variable. Ac- 
cordingly, we replace d m (r) by d(r, 9). For simplicity, we 
restrict ourselves to the case of broken time-reversal in- 
variance throughout the Brillouin zone. This corresponds 
to 7„ = 1 for all n. In order to enable a full quantitative 
comparison with Ref. jq] , we had to calibrate the dimcn- 
sionless conductance on basis of a common definition: In 
the present limit, and for x > 1, we find the following 
relation between the variance of the level velocities, the 
second moment of the K n {r) with respect to n, and the 
conductance, 



4tT 



' dr a (e y 

A9 



E Sr*»M = 9 -r- (is) 



afi n=~oo 



In Ref. H , the following expression is obtained for the 
flux-averaged correlation function 

(d(r,9)d(r + r,9 + 9)) Tl j-l = 



2tt 



9uc9 2 J A 



"dA 




exp 



2tt 



2tt 



,X + - 



ttA 



cos(rA) 



(19) 



An expression for K n (r) is reached by Fourier transfor- 
mation with respect to r and with respect to 9 over the 
Brillouin zone. The discrete spatial Fourier transform 
cannot be applied directly to Eq. ([l9]) since it violates 
periodicity in 9. Instead, exploiting the fact that this 
expression decays to zero for 9 — > 00, we approximate 
the Fourier sum over the Brillouin zone by a continu- 
ous Fourier transformation over the whole real axis. We 
obtain 



K n (r) = 



2g u 



■ {$n(x+) - <£>„(x-)} 



(20) 



with $„(x) = 7rnerf(n/2 v / a;) + 2y / 7nFexp(— n 2 /Ax) and 
x± = (2^)- 1 3uc |r±r 2 |. 

We have checked analytically and numerically (Fig. @) 
that both for r — > and r — > 00, Eq. ( P0[ ) coincides 
asymptotically with our results. Only in the vicinity of 
t = 1, there exist deviations (e.g., from Eq. ([17|) we find 
#o(l) = 1/V2ff^, while Eq. © gives K (l) = l/y/g^, 
cf. Fig. 0). The close agreement between the respective 
results provides further support for the matching proce- 
dure we used to connect the short- and long-time regimes. 
This agreement is not surprising, since it has recently 
been shown pl| -p3|l that the random-matrix approach 
underlying Eq. (M) also applies to systems where the 
disorder is of dynamical origin. 




FIG. 2. Comparison of the form factors K„(t) in the 
limit TV — ► 00, according to Eqs. (|10j) and (jL7|) (solid 
lines), to the corresponding functions as implied by the 
results of Ref. §), cf. Eq. ([§) (dashed). The graphs 
shown are for (top to bottom) n — 0, 16, 32, . . ., 256. 
The dimensionless conductance is g uc = 2007T, the same 
value as that underlying the data in Fig. Wa. 



VI. NUMERICAL RESULTS AND DISCUSSION 

In the remainder, we compare the theory sketched in 
the previous paragraphs to two prototypical yet quite dif- 
ferent models. Our first example is the kicked rotor on a 
torus ]24J] , defined by its Hamiltonian 



H(l,#;t) = 



(/-A) 2 



Va, h {&) 



OO 

E 



5(t-mT). (21) 



It is periodically time dependent, so that spectrum 
and eigenstates are adequately discussed in terms of 
quasienergies and Floquet states, respectively. The 
kicked rotor attains periodicity also in the angular- 
momentum variable I if the parameter r, an effective 
quantum of action, is chosen as r = Attp/ q, with p, q 



coprime (here we restrict them further to p = 1 and q 
odd). The unit cell then accomodates q quanta of angu- 
lar momentum. The number of quasienergy bands (for 
a typical sample, see Fig. ||) is also q. It is analogous to 
the total number of bands Ar introduced above. 

The quasienergies are obtained by diagonalizing the 
symmetry-projected Floquct operator at Bloch phase 9 m 



(l'\U m \l)=exp(-2n£(l-\) 2 



1 9_1 

7X"- 1 



V Q , fe ([0 m +27rn]/qr) Ml-l')(6 m +2Trn) / q 



(22) 



n=0 



For the kicked rotor on a torus, the number of unit cells 
is simply determined by the number N of Bloch phases 
where Eq. d23) is evaluated. This additional parameter is 
independent of q, the number of bands. The total num- 
ber of levels in the spectrum is therefore Nq. 




FIG. 3. Quasienergy bands pertaining to chaotic 
Bloch states of the quantum kicked rotor, for k = 300 
and t — 4-7T/75, corresponding to a periodic potential 
with a unit cell that accomodates 75 Bloch states. 

Besides being periodic in I and V , the Floquet ma- 
trix ( p2|) allows for several twofold symmetries P4|-j26|. 
In order to break them in a controlled manner without 
significantly altering the classical diffusion, an angular- 
momentum shift by A has been introduced, and the po- 
tential has been chosen as 129] 



V a A§) 



IT 1 7T 

cos(o;— ) cos $ H — sin(a— ) sin 2i3 



(23) 



For a = A = m = 0, which classically corresponds to 
the common standard map, the kicked rotor possesses 
two independent twofold symmetries, a unitary one, the 
parity P: I — > —I, $ — ► 2tt — i9, t — > t, and an antiunitary 
one, time reversal T: I — > — I, ?9 — > ?9, t — >• —t. P invari- 
ance is broken for m ^ 0, L/2, i.e., off the band center 



and edges. Therefore, for a = A = 0, the quasienergy 
statistics corresponds to the superposition of two inde- 
pendent COE's at m — 0, ±L/2, and to a single COE 
elsewhere. Choosing a ^ breaks P also at the symme- 
try points and at the same time, lifts the band symme- 
try so that COE statistics becomes valid throughout the 
Brillouin zone. If A takes a non-integer value, T invari- 
ance, which amounts to (—l'\U m \ — I) = (l'\U m \l), is 
broken. In this case, the quasienergy statistics is that of 
the CUE everywhere, or that of two superposed CUE's 
at the symmetry points if a — 0. 
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FIG. 4. (a) Form factors K u (t) for the kicked ro- 
tor on a torus at (top to bottom) n = 0, 16, 32, . . ., 
256 (N = 512), compared to the theory (heavy lines). 
The parameters determining <; uc are, in the notation of 
Ref. Ill, r = 4-7T/225 (dimensionless quantum of action, 
not to be confounded with the scaled time used in this 
letter) and k — 300, giving g uc = 200-7T. All twofold an- 
tiunitary symmetries are broken, (b) is analogous to (a) 
but for a billiard chain composed of N — 128 unit cells 
(inset). The winding numbers shown are n = 0, 2, 4, 8, 
16, 32, 64. The only free parameter of the theory, g uc , was 
determined by an independent classical simulation Jl8| . 

Sufficiently large spectral data sets have been gener- 
ated by varying k over small intervals so that the classical 
diffusion constant D rs fc 2 /2 is not changed appreciably. 
The parameters N, q, and k were chosen such that the 
scale of localization due to disorder within the unit cell 
[p|,p4| by far exceeds the cell size. In Fig. Ha, we show, for 



both twofold symmetries simultaneously broken (a = 0.2, 
A = 0.5), the set of K n (r) as a function of r at selected 
equidistant values of n, and compare with our theory, 
Eqs. ( |10|) and (|16|). The three time regimes and the 
corresponding power-law time dependences of Kq(t) can 
clearly be distinguished. Comparing with the nonlincar- 
sigma-model results [pi (Eq. ( |20| ) and the dashed lines 
in Fig. @), we find that near r = 1, they reproduce the 
data even better, including an unexpected feature like 
the sharp peak at n = 0, r = 1. The saturation regime 
r ~ N^n/g uc , however, is not contained in Eq. (p0|) . 

Our second example is a billiard chain composed of 
unit cells as shown in the inset in Fig. [|b. The wave 
functions satisfy the Hclmholtz equation augmented with 
Dirichlet boundary conditions on the channel walls and 
periodic boundary conditions along the channel after N 
unit cells. Provided that trajectories traversing the unit 
cell without hitting the obstacles are excluded by an ap- 
propriate geometry, the classical dynamics is diffusive on 
scales larger than the size of the unit cell. The energy 
bands can be found using the scattering approach to bil- 
liard quantization, see Ref. Jig] for details of the method. 
The shape of the unit cell is chosen such that there are no 
other unitary symmetries besides the discrete translation 
invariance. However, due to time-reversal invariance, the 
bands are symmetric with respect to 6 = and it, such 
that 70 = 7at/2 = 2 and else 7„ = 1. Fig. [|b, which 
is analogous to Fig. Ha, clearly shows the correspond- 
ing enhancements of Kq(t) and K^nir) throughout the 
three regimes. Correspondingly, the spectral statistics 
for fixed Bloch number follows the COE at the center 
and the edge of the Brillouin zone and approaches CUE 
in between p8[ . 

Our theoretical results reproduce remarkably well the 
numerical data obtained for the two models. This lends 
support to the approximations underlying the theory and 
corroborates its conclusions. The large- AT limit is of 
particular interest, because it provides a very sensitive 
tool for studying the transition to Anderson localiza- 
tion when the periodicity is disrupted by disorder. In 
the latter case, the level statistics approaches the Pois- 
sonian limit, and the corresponding spectral form factor 
approaches the constant value 1 in the long-time limit. 
This is in a marked difference from the periodic case, 
where Kq(t) « 1/t for r > l/N. The transition between 
these two rather distinct asymptotics occurs as a smooth 
function of the degree of disorder. In order to understand 
this transition in semiclassical terms [0,0, one may dis- 
cuss the effect of weak disorder through the small ran- 
dom differences between the phases assigned to periodic 
orbits which were related by symmetry in the absence 
of disorder. This approach is similar to the semiclassical 
discussion of the response of the weak-localization peak 
to a weak magnetic field pj| . A detailed analysis of this 
scenario will be left for further study. 



ACKNOWLEDGEMENTS 



The research reported in this work was supported by 
grants from the Minerva Center for Nonlinear Physics. 
Three of us (TD, BM, HS) would like to thank the Weiz- 
mann Institute of Science, Rehovot, and US would like 
to thank the Max Planck Institute for Physics of Com- 
plex Systems, Dresden, for the kind hospitality enjoyed 
during several visits. 



[1] T. Dittrich and U. Smilansky, Nonlinearity 4, 85 (1991); 
T. Dittrich, E. Doron, and U. Smilansky, J. Phys. A 27, 
79 (1994); T. Dittrich, Phys. Rep. 271, 267 (1996). 
[2] N. Argaman, Y. Imry, and U. Smilansky, Phys. Rev. B 

47, 4440 (1993). 
[3] J. H. Hannay and A. M. Ozorio de Almeida, J. Phys. A 

17, 3429 (1984). 
[4] M. V. Berry, Proc. R. Soc. A 400, 229 (1985). 
[5] B. D. Simons and B. L. Altshuler, Phys. Rev. Lett. 70, 

4063 (1993). 
[6] N. Taniguchi and B. L. Altshuler, Phys. Rev. Lett. 71, 

4031 (1993). 
[7] T. Geisel, R. Ketzmerick, and G. Petschel, Phys. Rev. 

Lett. 66, 1651 (1991); 67, 3635 (1991); 69, 695 (1992). 
[8] E. R. Mucciolo, R. B. Capaz, B. L. Altshuler, and 

J. D. Joannopoulos, Phys. Rev. B 50, 8245 (1994). 
[9] H. Silberbauer, P. Rotter, U. Rossler, and M. Suhrke, 
Europhys. Lett. 31, 393 (1995). 
[10] G. A. Luna-Acosta, K. Na, L. E. Reichl, and A. Krokhin, 

Phys. Rev. E 53, 3271 (1996). 
[11] A. R. Kolovsky, S. Miyazaki, and R. Graham, Phys. Rev. 
E 49, 70 (1994); A. R. Kolovsky and S. Miyazaki, Phys. 
Rev. E 50, 910 (1994); A. R. Kolovsky, Phys. Rev. E 56, 
2261 (1997). 
[12] H.-J. Korsch, B. Mirbach, and B. Schellhaafi, J. Phys. A 

30, 1659 (1997). 
[13] J. M. Robbins, Phys. Rev. A 44, 2128 (1989). 
[14] S. C. Creagh and N. D. Whelan, Phys. Rev. Lett. 77, 

4975 (1996). 
[15] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics, Interdisciplinary Applied Mathematics, vol. 1, 
Springer (New York, 1990). 
[16] P. Leboeuf and A. Mouchet, Phys. Rev. Lett. 73, 1360 

(1994). 
[17] R. Scharf and B. Sundaram, Phys. Rev. Lett. 76, 4907 

(1996). 
[18] T. Dittrich, B. Mehlig, H. Schanz, and U. Smilansky, 

Chaos, Solitons & Fractals 8, 1205 (1997). 
[19] J. R. Stine and R. A. Marcus, J. Chem. Phys. 59, 5145 

(1973). 
[20] N. W. Ashcroft and N. D. Mermin, Solid State Physics, 

Holt-Saunders (Philadelphia, 1976). 
[21] B. A. Muzykantskii and D. E. Khmelnitskii, JETP Lett. 
62, 76 (1995). 



[22] O. Agam, B. L. Altshuler, and A. A. Andreev, Phys. Rev. [25] R. Bliimel and U. Smilansky, Phys. Rev. Lett. 69, 217 

Lett. 75, 4389 (1995). (1992); M. Thaha, R. Bliimel, and U. Smilansky, Phys. 

[23] E. B. Bogomolny and J. P. Keating, Phys. Rev. Lett. 77, Rev. E 48, 1764 (1993). 

1472 (1996). [26] G. Casati, R. Graham, I. Guarneri, and F. M. Izrailev, 

[24] F. M. Izrailev, Phys. Rep. 196, 299 (1990). Phys. Lett. A 190, 159 (1994). 



